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Gravitational lensing of distant objects caused by gravitational tidal forces from inho- 
mogeneities in the universe is weak in most cases, but it is noticed that it gives a great 
\^ ■ deal of information about the universe, especially regarding the distribution of dark mat- 

' ter. The statistical values of optical quantities such as convergence, amplification and shear 

have been derived by many people using various approaches, which include the linear per- 
turbational treatment in the weak limit and the nonlinear treatment considering small-scale 
^ ■ matter distribution. 

' In this review paper we compare the following three main approaches: (a) the approach 

\ in the multi-lens-plane theory; (b) the approach due to the direct integration method; and 

^f) , (c) the perturbational approach. 

' In the former two approaches inhomogeneous matter distributions are produced in the 

' CDM model using A'^-body simulations (the P'^M code and the tree-code, respectively). In 

^\ , (c) the power spectrum corresponding to the CDM model is used for the large-scale matter 

^\ . distribution. 

Q-i| §1. Introduction 

o . 

^ , The propagation of light from distant objects Hke galaxies and QSOs is deflected 

^ I by the gravitational tidal forces caused by inhomogeneous matter distribution be- 

tween the objects and us. The so-called lens effect creates conspicuous images such 
as multiple QSOs and arcs in clusters of galaxies, owing to special positional relations 
^ ' among sources, lens objects and the observer, but in most cases it causes small de- 

■ formations and amplification of optical images, which result from the superposition 

of deflections from many lens objects on cosmological scales. This so-called weak 
lensing gives us valuable information on the structure and evolution of the universe, 
especially regarding the distribution of dark matter, not only on large scales, but 
also on small scales around galaxies. 

The statistical behavior of optical quantites such as convergence, amplification 
and shear due to weak leasing has been studied by many people since the pioneer- 
ig papers by Gunn, Weinberg cP and Blandford and Jaroszyhski.cl) For the 
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derivations of these quantities there have been various approaches which consist of 
the multi-lens-plane method, direct integration methods solving the null-geodesic 



E-mail address: tomita@yukawa.kyoto-u.ac.jp 

E-mail address: premadi@astr.tohoku.ac.jp 

E-mail address: nakamura@utaphp2.phys.s.u-tokyo. ac.jp 



86 



K. Tomita, P. Premadi and T. T. Nakamura 



equation, optical scalar equations and the equation of geodesic deviation, and per- 
turbative methods. 

1.1. Multi-lens-plane method 

This is one of the various methods numerically simulating light propagation in 
inhomogenous model universes, whicli^were derived and developed by BlaridforiLand 
Narayan,!^ Schneider and Weiss, EP"@) Jaroszyhski et al.,EP Jaroszyhski.ty^ "IHI' In 
this method we first build a finite number of planes normal to light rays between 
an observer and a source, and the three-dimensional distribution of galaxies and 
dark matter as lens objects is replaced by their two-dimensional distribution in the 
planes onto which the matter in each interval is projected. The deflection of light 
rays resulting from each lens plane is calculated using geometrical optics^ and the 
statistical averages of optical quantities for many-ray bundles are derived. 

To use this method we must produce the inhomogeneous distribution of lens 
objects in each lens plane. In earlv works (Schneider and Weiss, eP Paczyhski and 
Wambsganss,llHP Lee and Paczyhski. the random distribution of lens objects with 
an equal mass was assumed, and so the realistic large-scale structure of the matter 
distribution was neglected. Jaroszyhski et al.i^ constructed an inhomogeneous model 
using an A^-body simulation (the PM code), but in this model large-scale structure 
was consideped, j^hile small-scale structure corresponding to galaxies was neglected. 
Jaroszyhski 113 "EJ) developed other numerical models of matter distribution using 
the Zeldovich approximation, located galaxies (as lens objects) correponding to the 
matter density, and assumed for them the mass spectrum due to the Schechter 
luminosity function and their morphological types randomly. 

Wambganss et al.0 "@ generated models for the evolution of large-scale struc- 
tures by A^-body simulation in the CDM models, investigated their influence on light 
propagation and statistical results such as the frequency and separation angles of 
multiple QSOs and the image deformation of high-redshift objects. In their works, 
galaxies as lens objects were not taken into account. 

Recently Premadi et al.llB' have improved the work of Jaroszyhski et al. by 
adopting an N-hody simulation (P'^M code) and considering the spatial distribution 
of galaxies with the mass spectrum and morphological types (due to the morpholog- 
ical type-density relation). They adopted Jaroszyhski's assumption that the main 
lens objects are galaxies and the background matter outside galaxies has smooth 
distribution with radii ~ l/i~^ Mpc. In §2 the recent result of studies by Premadi 
et al.x3 is given. 

1.2. Direct integration methods 

There are three types ofdirect integration methods in which the optical scalar 
equation (derived by Sachs EJ)), the equation of geodesic deviation and the null- 
geodesic equation are solved. Kantowski's approach belongs to the first type. He 
derived the averaged optical scalar equation and solved it in an inhomogeneous model 
universe with swiss cheese structure. Dyer and Roeder© studied the behavior 
of the amplification applying this aijproach to the case of a nonzero cosmological 
constant. Recently, Kantowski et al.B2P have analyzed the lens effect in th£_pbserved 
redshift-magnitude relation for type la supernovae. Watanabe and TomitaclP studied 
the behavior of shear in the ray bundles in various models in which the random 
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distribution of particles simulating galaxies and clusters of galaxies was assumed for 
simplicity. As an example of the second type of direct integration method (treating 
the equation of geodesic deviation) we cite a recent work of Holz and Wald,E2P 
who studied the behavior of amplification in an inhomogeneous model with a swiss 
cheese-like structure. In their case not all systems consisting of a central galaxy and 
the surrounding empty region have the compensated mass distribution, in contrast 
to Kantowski's swiss cheese system. They derived the probability distribution of 
amplification, which was found to be non-Gaussian in accord with Dyer-OattesE^' 
and examined the significance of the amplification on the observation of type la 
super novae (Holz§) ) . A drawback of this approach may be that it treats only ray 
bundles with infinitesimal separation angles. 

The third tvne of direct integration method is described in the works of Tomita 
and Watanabe.c3) "Ej) They treated the multiple deflection of light rays directly 
by solving the null-geodesic equation and studied the statistical averages of image 
deformation in inhomogeneous universe models. The evolution of models was derived 
under the periodic condition using Aarseth's individual particle method for A^-body 
simulations, in which particles are regarded as galaxies or clusters of galaxies with 
the initial power spectrum n = 0. The particles (with the same mass and particle 
number = 1331) are placed and move in a periodic box with comoving length 
~ 50 Mpc, and light rays also propagate in the box under gravitational forces from 
particles in the box. In Tomita and Watanabe,© the influence of lensing on the CMB 
anisotropy was analysed^jida negativejiesult was obtained, contrary to expectations 
(Kashlinsky,KlP Tomita,E§) "Ea) Sasaki. eS)) In the latter paper, liB' examples of strong 
image drfomiations and their statistics were studied. Using the same approach, 
Tomita e1|) 'c3 has recently studied the statistical behavior of optical scalars in more 
realistic inhomogeneous models with the CDM spectrum (n = 1). All particles are 
assumed to have the same galactic mass, while their radii depend on the model used. 
In the first two of these papers, "EHP the lens effect was analyzed in flat and open 
models, respectively, and in the third paper, Ej) the distribution of angular diameter 
distances in these models was derived. 

As for the number density of main lens objects, there is an important differ- 
ence between the work of Premadi and that of Tomita. In the former, adopting 
Jaroszyiiski's assumption mentioned above, the density of main lens objects (galax- 
ies) is found to be f2g ~ 0.02. In Tomita's works, on the other hand, the density 
{f^L) of main lens objects (called compact lens objects) is i?^ = 0.2 for three models 
with Qq = 1 and 0.2 and i?^, = 0.4 for a model with Qq = 0.4, which are much 
larger than Qg. This difference between these two works comes from the differ- 
ent estimation for the role of non- galactic clouds as lenses. At present it is not clear 
whether Jaroszyhski's assumption is realistic or not, and so the matter outside galax- 
ies may not be so smooth and may behave as discrete clouds with radii ~ 100 
kpc, consisting of dark matter and the baryon fluid. In these clouds a baryon fiuid 
has experienced no dissipated collapse and has less central concentration than in 
galaxies, but as weak lens objects they may play a role similar to galaxies. 

An outline of Tomita's recent works is given in §3, and the results of additional 
analyses are also given in the case (satisfying Jaroszyhski's assumption) that the 
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density of lens objects is Ql = and the remaining objects are assumed to be 
smoother, by constraining their radii so as to be 500/i~^kpc. By comparing the two 
results in the case of Ql = 0.2 and the case of Ql = ^g-, we find how lensing on the 
small scale is dominated by non-galactic lens objects, if they exist. 
1.3. Perturbative methods 

The method for solving optical scalar equations perturbatively with respect to a 
small perturbed expansion and shear has hem introduced by Gunn& and developed 
by Babul and LeeCS and Blandford et al.Ej). Gravitational influences from large- 
scale structures are involved in the terms giving the Ricci and Weyl focusing, through 
the density power spectrum or the numerically simulated matter distribution. The 
angular correlation functions of optical scalars were derived in connection with the 
two-point density correlatkm function, and higher correlations also were studied by 
Villumsen,E9 BernardeauE3) and Nakamura. H 

The perturbative formulation for the image deformations has been derived inde- 
pendently from the null-geodesic equation (Linder,L3 Martinez et al.,E2P Martinez 
and SanzeiP). This approach also has been applied to studies concerning weak im- 
age deformation of high-redshift objects and thesmoothing of the CMB anisotropy 
(Cayon et al.,© Seljak© and Jain and SeljakcS). In §4, a recent result obtained 
using the power spectrum approach by Nakamura is given. 

§2. Lensing in the mult i- lens- plane method 

2.1. Methodology 

The inhomogeneities that perturb a light bundle as it travels from a source to an 
observer are of sizes much smaller than the distance that the light traverses. This has 
motivated a number of authors to use the multi-plane gravitational lensing approach 
to study how light propagates in inhomoeeneous universes (e.g. Schneider and Weiss, 
Jaroszyiiski et al.,^' Jaroszyhski,0* "0^ Wambsganss, Cen and Ostrikertl^). 
In this approach one first idealizes the inhomogeneities as being distributed on a 
series of thin lens planes which are arranged perpendicular to the line of sight. Then 
one assumes that lensing only occurs in each of those planes. This way one can 
analyze the lensing properties of each plane separately and let the light beam carry 
the lensing effect of each plane while propagating from one plane to the next. For 
this paper to be self-contained we give a brief review of multi-plane lensing theory, 
following closely the description and notation of Schneider et al.W 

2.2. Multi-plane gravitational lensing 

Consider N lens planes located at redshifts Zi, with i = l,N, and ordered such 
that Zi < Zj for i < j. Figure 1 displays an example with N = 2. All angles are 
exaggerated for clarity. Each lens plane is characterized by its respective surface 
mass density where is the impact vector of the ray in the i-th lens plane. 

Let cuii^i) denote the deflection angle the light ray experiences in the i-th plane at 
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a position ^j. From this geometry, we can derive the lens equation, 

n ^ 

^ = :^^i-Ea5a.(^j, (2-1) 
^ 1=1 

where rj is the source position vector (in the source plane), is the impact vector in 
the i-th plane, Dj is the angular diameter distance between the j-th plane and the 
observer, and Dij is the angular diameter distance between the i-th and j-th planes, 
with S = N + 1 identifying the source plane. Knowing the impact vector in 
the image plane, the impact vector in subsequent planes can be obtained recursively 
using 

D 

^i = 77«i-EAA(^.)- (2-2) 

1 1=1 

The deflection angle is related to the surface density by 

"^(0 = ^ // , (2-3) 

where G is the gravitational constant, c is the speed of light, and the integral extends 
over the lens plane. We can rewrite this expression conveniently as 

"^(^J = V^.(^,), (2-4) 

= ^JJ ^^(^0 In 1^. - ^'M'C' • (2-5) 

It is useful to rewrite these equations in a dimensionless form. Define for each lens 
plane a critical surface density as 

c^Ds .... 

and introduce the following dimensionless quantities : 

a;» = |f, l<i<iV+l; (2-7) 

K.[xi) = —, l<i<N. (2-8) 
Equations (2-4) and (2-5) reduce to 

^3 = =^1 - E PijO^ii^i) , (2-9) 

i=l 

ai{x^) = ViPi{xi), (2-10) 
'ilji{xi) = — JJ Ki{x') ln\xi - x'\(fx' , (2-11) 

where 
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Fig. 1. Multi-lens-planes. 

and the gradient is now taken relative to Xi. We can invert equation (2-11) and 
obtain 

V'^'^|J^ = 2Ki. (2-13) 

To compute the scaled position y = xs oi the source in the source plane, we simply 
set j = N + 1. Equation ( p-12| ) gives I3is = 1, and Eq. (2-9) becomes 

N 

y = XN+i = xi-^OLi{x.i). (2-14) 

i=l 
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This ray-tracing equation is a mapping from the image plane (i = 1) onto the source 
plane {i = N + 1). 

The effect of each lens plane on the evolution of the beam is described by the 
Jacobian matrix 

^.(-.)=r^ (2.15) 



-tpi,21 1 - -0^,22 . 

where the commas denote differentiation with respect to the components of Xi. Since 
■04,12 = '0j,2ii and Eq. ( |2-13| ) gives tpi^u + V'i,22 = 2«;j, we can rewrite Eq. (2-15) as 

. f 1 — Ki — Sn —Si2 \ 1«^ 



where 



We now define 



511 = ^(V'i.ii-te), (2-17) 

512 = i^^^2 = A,21 . (2-18) 



S, = {Sl + S!,y/\ (2-19) 



The determinant and trace of Aj can be expressed entirely in terms of Ki and Si, as 
follows: 

det A, = (1 - Kif - Sf , (2-20) 
trAi = 2(1 - Ki) . (2-21) 

The quantities fii = l/(det Aj), 1 — Ki, and Si are called the magnification, conver- 
gence (or Ricci focusing), and shear, respectively. 

To compute the cumulative effect of all the lens planes, we consider the Jacobian 
matrix of the mapping given by Eq. ( 2T4 ), 

Bix) = ^ = I-Y^pl = I-Y: U.B, , (2-22) 
where / is the 2x2 identity matrix, and Ui and Bi are defined by 

U^ = ^, (2-23) 

OXi 

B. = ^. (2-24) 
ox I 

After substituting Eq. (2-10) into Eq. (2-23), we obtain 

t/. = /-A,= f^^-i^ ^^'12 V (2-25) 

\^i,2i ^i,22 y ' ^ ^ 

where Aj is given by Eq. ( 2-15| ). Hence, Ui describes the effect the i-th plane would 
have on the beam if all the other planes were absent, and Eq. ( 2-22| ) simply combines 
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the effect of all the planes. To compute the matrices -Bj, we differentiate Eq. (2-9) 
and get 

i-i 

Bj = I-Y, AjC/^Si . (2-26) 

i=l 

Since Bi = I, we can use Eq. ( |2-26D to compute all matricies Bi by recurrence. 

2.3. The numerical algorithm 

We use the P^M algorithm (Hockney and Eastwood©) for the A^-body simu- 
lations of the large scale structure (LSS) of the universe. The calculations evolve a 
system of gravitationally interacting particles in a cubic volume with triply periodic 
boundary conditions, comoving with the Hubble flow. The forces on the particles 
are computed by solving the Poisson equation on a cubic grid using a Fast Fourier 
Transform method. 

2.3.1. Building the lens planes 

To implement the multi-plane lens method, we divide the space between the 
source and the observer into a chain of cubic boxes of equal comoving size, Lbox- We 
first need to determine the redshifts of the interfaces between these cubic boxes. Let 
us assume that the photons that reach the observer at present entered a particular 
box at time t' and redshift z' and exited that box at time t and redshift z. The 
redshifts z' and z are related by (Premadi et al, 



L 



box 



t' 



l + z{t) cdt. (2-27) 



Using this equation, with the appropriate relation for z{t), we can find the redshifts 
of the interfaces. The front side of the box closest to the the observer is, by definition, 
at z = 0. Plugging this value into Eq. ( p-27 ) gives us the redshift z' of the back side 



of the box, which is also the redshift z of the front side of the next box. Then, by 



using Eq. ( 2-27 ) recursively, we can compute the redshifts of all the interfaces. 

As the structures inside the box might evolve while the light beam travels across 
it, we decide the plane onto which we project the mass distribution to be at redshift 
where the density contrast is equal to the time-averaged density contrast of the box. 

2.3.2. The galaxy distribution 

We consider the LSS at present {z = 0) resulting from the P^M simulations, and 
we design an empirical Monte-Carlo method for locating galaxies in the computa- 
tional volume, based on the constraints that (1) galaxies should be predominantly 
located in the densest regions, and (2) the resulting distribution of galaxies should 
resemble the observed distribution on the sky. Our method is the following. We 
divide the present computational volume into 128^ cubic cells of size IMpc^, and 
compute the matter density p at the center of each cell using the same mass assign- 
ment scheme as in the P'^M code. We then choose a particular density threshold pf 
We locate N galaxies in each cell, where is given by 

N = mi(^j-^. (2-28) 
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The actual location of each galaxy is chosen to be the center of the cell, plus a 
random offset that is of the order of the cell size. This eliminates any spurious effect 
introduced by the use of a grid. We then experiment with various values of the density 
threshold pt until the total number of galaxies comes out to be of order 40000. This 
gives a number density of ~ 0.02^alKxies /Mpc^. This method bears some similarities 
with that used by Jaroszyiiski. t3) 't3 Tests showed that the observed galaxy 2-point 
correlation function is fairly well reproduced (Martel et al.E^P). 

The morphological type of each galaxy is determined by using a combination of 
the known relations between the present distribution of morphological types and the 
surface density of galaxies along with a Monte-Carlo method. 

The locations of each galaxy at higher redshifts are determined by combining 
the distribution of the galaxies and that of the particles from the P^M simulation, 
i.e., by following the position of the nearest particle. 

We adopt the galaxy models described in Jaroszyhski. The projected 

surface density of each galaxy is given by 

a{r) = I 4(7(^2 + ^2)1/2 ' ^ < ^"1^^ ' (2-29) 

I , r > Tmax , 

where r is the projected distance from the center. The parameters , ?^max and 
V are the core radius, maximum radius, and the rotational velocity, respectively. 
These parameters are functions of the luminosity and morphological types of each 
galaxy. We assume that the present galaxy luminosities are given by the Schechter 
luminosity function, 

n{L)dL = ^ (IlY e~L/^'dL, (2-30) 

where n[L) is the number density of galaxies per unit luminosity, a = —1.10, n* = 
0.0156/i^Mpc~^, and Lg* = 1.3 x lO^°/i~^L0, where Lb is the luminosity in the B 
band (Efstathiou et al.O). For each galaxy we generate a luminosity L according to 
this distribution, and we combine it with the galaxy morphological type to determine 
the values of the parameters f , Tc and rmax for that galaxy. 

2.4. The experiments 

2.4.1. The cosmological models 

We are currently appyling the algorithm described above to perform a large cos- 
mological parameter survey. We focus on the CDM model whose density fluctuation 
power spectrum is described in detail in Bunn and White. This power spectrum 
is characterized by 6 independent parameters : (1) the density parameter Qq; (2) 
the contribution Qb of the baryonic matter to the density parameter; (3) the cosmo- 
logical constant Aq; (4) the Hubble constant Hq; (5) the temperature Tcmb of the 
cosmic microwave background radiation; and (6) the tilt n of the power spectrum. 
We initially set Tcmb = 2.7K and Qb = 0.015/i~^, thus reducing the dimensionality 
of the parameter space to 4. The normalization of the power spectrum is often de- 
scribed in terms of the rms density fluctuation erg at a scale of 8/i~^Mpc, which is a 
function of the 6 aforementioned parameters. We invert this relation, treating erg as 
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an independent parameter, and the tilt n as the dependent parameter. Presently we 
have 43 models with various (but restricted) values of the above first 4 parameters. 

2.4.2. The ray shooting 

In each model we propagate beams of 341^ light rays backward to redshift Zg ~ 3. 
The beam has a size of 21.9 arcseconds, and the spacing between rays is 0.064 
arcseconds. To analyze the results of the experiments, we lay down grid on the 
source plane, divided into 31 x 31 cells. Cells which have more than the average 
number of rays (34lV3l2=121) indicate sources which are magnified. The fraction 
of such cells gives us the magnification probability. By taking the rays located in 
such cells, and tracing them back to the image plane, we can study the shape of the 
images, and in particular, compute the fraction of the magnified sources that have 
double images. Figure 2 displays a few examples of images of a circular source. 

2.5. Preliminary results 

We have so far conducted between 30 and 60 ray shooting rays on all of the 
models, but results of only 35 models have been analyzed. In a previous work 
(Premadi et al. 1998) we tested the methodology we have described here by studying 
the statistics of magnification and shear with respect to the lens redshift for a given 
source redshift. Among other things, we have found that the lensing effect is most 
prominent at intermediate rcdshifts although the structures are more evolved at lower 
redshifts. The redshift where most lensing occurs depends upon the cosmological 
models. We also show that the magnification is dominated by convergence, with 
shear contributing less than one part in 10^. 

In the current work we study the magnification probability, Pm, the probability 
of double-image events given a lower limit of magnification, P2, and the distribution 
of the separation angle in the double image events. 

In Fig. 3, despite the still insufficient statistics, for the majority of models, we 
observe a tendency for P„j to decrease with increasing erg, and for the cases in which 
we do not observe this tendency, it is not ruled out because of the large error bars. 
Having an anticorrelation between P^ and a% is seemingly counter-intuitive, since erg 
measures the amplitude of the density fluctuation, which are responsible for lensing. 
The explanation is that the magnification is caused primarily by the matter located 
near the beam, whereas matter located far from the beam are responsible for the 
shear. In cosmological models like CDM, structure formation proceeds hierarchically. 
Small structures form first, then merge to form larger structures. A larger cjg implies 
that this hierarchical merging process is more advanced. This means that the clusters 
are more massive, and possibly also denser, but there are fewer of them, thereby 
enlarging voids between them. The light beam is therefore less likely to hit or 
pass near any cluster, resulting in a smaller magnification probability for larger cjg. 
However, if the beam does hit a cluster, we would expect the effect to be stronger in 
models with large erg, as the clusters are more massive. 

In Fig. 4, the plot of P2 versus P^ indeed shows an anticorrelation for several 
models. Another interesting result is that at fixed Qq, Aq, and erg, the magnification 
probability is insensitive to the value of Hq. This results from the combination of 
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(a) 



magnification and shear 



(c) 













Einstein ring 



(b) 

















doubie image 



(d) 



tripie image 



Fig. 2. Various images of circular sources. 



two competing effects. On one hand, the number of lens planes between the source 
and the observer decreases with increasing Hq. On the other hand, the total mass in 
each plane is proportional to the critical density, and thus increases with Hq. Hence, 
models with larger Hq have a smaller number of more massive lens-planes. 

The histograms of the separation angle (Fig. 5) shows the presence of double 
peaks in most cases. This seems to be more prominent in the models with larger cxg 
when the other cosmological parameters are held fixed. The high density clusters 
might contain more elliptical galaxies which have smaller cores than the field galaxies. 
A more compact core is known to result in a larger separation angle. Intriguingly, 
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Fig. 3. Magnification probability Pm of magnification 1.2 and larger as a function of the rms den- 
sity fluctuation a%. The numbers in parentheses indicate the values of (i7o,Ao). Open circles 
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Fig. 4. Multi- image probability P2 as a function of magnification probability Pm- The numbers in 
parentheses indicate the values of (J7o, Ao). Open circles correspond to Ho = 65 and circles to 
Ho = 75. 



.12 
.1 
.08 

qT .06 
.04 
.02 


I 

.1 
.08 
.06 

CM 

Q. 

.04 

.02 


.1 

.08 

.06 

CM 

Q. 

.04 
.02 


I 

.1 

.08 

.06 

CM 

CL 

.04 
.02 




T — \ — I — I — I — \ — I — I — I — \ — I — I — I — \ — I — I — \ — I — r 



O 



O 



- (0.2, 0.0) 

J \ I I I \ I I I \ I I I \ I I \ I L 



.15 



.05 



T — I — I — \ — I — I — \ — I — r 
O ' 
CD 



(0.2, 0.8) 

J \ I I I I I I \ I I \ I L 



.02 .04 .06 .08 .1 

I I I I I I I I I I I I I I I I I I I I I I I 

; • 

o 



o 



.05 



.1 



.15 



(0.5. 0.0) 



I I I I I I I 



I I I I I I I I I I 



.02 .04 .06 .08 .1 .12 



_ 1 1 1 1 1 


1 1 1 1 

' : 




• — 




• 




— 


~ (0.7, 0.0) 




" 1 1 1 1 1 


1 1 1 1 1 " 



.1 

.08 

.06 

.04 

.02 


.1 

.08 

.06 

.04 

.02 



I I I I I I I I I I I I I J I I I I 



O 



(0.5, 0.5) 



J I I \ I I I \ I I \ I I I \ I I L 



.02 .04 .06 .08 .1 



T — I — I — \ — I — I — I — \ — I — I — ^ — I — I — I — ^ — I — I — r 
O _ 



O 



(0.7, 0.3) 



.02 



.04 



1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 






• 


° ° J 




• 


~ (1.0, 0.0) 




~ 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 ~ 



I I I I I I I I I I I I I I I I I I 
06 "0 .02 .04 .06 .08 .1 



.02 .04 .06 .08 .1 



98 



K. Tomita, P. Premadi and T. T. Nakamura 



Fig. 5. The distribution of the image separation in arc seconds. The counts are not normahzed to 
the numbers of runs which are not the same for all models. Top row displays results for models 
with i7o = 1.0, Ao = 0.0. The middle row for models with Qo = 0.2, Ao = 0.0, and the bottom 
row for models with f2o = 0.2, Ao = 0.8. 
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models with smaller Hq seem to produce longer high separation tails. The presence 
of Aq also seems to stretch the histogram to a higher separation tail, although by a 
smaller amount than that caused by high erg. In addition to this, it also increases 
the counts at mid-separation. This might indicate that structures in Aq models have 
large sizes which are responsible for midsize separation angles. 

Even at this preliminary stage, the results show the different tendencies among 
various cosmological models. With better statistics, we are hopeful that the dif- 
ferences can be sharpened, and that these differences will eventually enable us to 
make meaningful comparisons with observational results, thereby discriminating the 
cosmological models. 

§3. Cosmological lensing in the direct integration method 

In this section we study the behavior of light ray bundles in inhomogeneous 
model universes, solving directly the null geodesic equation in the (cosmological) 
Newtonian approximation at the stage from an epoch ti (with redshift zi = 5) to the 
present epoch to- Light rays are deflected by lens objects such as galaxies and non- 
galactic clouds. Here, these clouds are mass concentrations which are unluminous 
but have masses comparable with a standard galactic mass. It is not clear at present 
what average mass and size and what number density such clouds have, compared 
with those of galaxies. In some previous papers (by Tomitail)-!!)), we assumed 
that non-galactic clouds are dominant and have a standard galactic mass and lensing 
strength similar to that of galaxies. In this paper we consider another case, in which 
only galaxies are dominant lenses and all non-galactic clouds are very weak lenses, 
and compare the results concerning optical deformations in this case with those in 
the previous case. In the future observations about the shear deformation at small 
angles will provide the information for the structure and lensing strength of non- 
galactic clouds. 

3.1. Model universes and the ray shooting 

The background model universes have the line-element 

ds'^ = -(1 + 2ip/c'^)c^dt^ + (1 - 2ip/c^)a\t){dxf/[l + K^{xf]'^, (3-1) 

where K is the signature of spatial curvature (±1, 0). The normalized scale factor 
S = a(t)/a(to) satisfies 



/dS_y _ 1 
[d^J ~ S 



Qo - {^0 + Ao - l)S + XoS^ 



(3-2) 



where r = Hot and ao(= a(to)) is specified by the relation {cHq ^/ao)^ = 1 — ■i^o ~ Aq. 
The gravitational potential ip is described by the Poisson equation 



a-'A^ = [1 - -{x)'f 



d'^ip ^ dtp 



dx^ [l + i^i(cc)2]3 5x^ 
A7rGpB[p{x)/pB-l], (3-3) 
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(3-4) 



where Hq = 100/iMpc^ km . In our treatment the inhomogeneities are locahy 
periodic in the sense that the physical situation at x is the same as that at a; + /n, 
where the components of n(= (n^,n^,n^)) are integers. In an arbitrary periodic 
box with coordinate volume Z^, there are particles with the same mass m. It is 
assumed that the force at an arbitrary point is the sum of forces from particles in 
the box whose center is the point in question, and that the forces from outside the 
box can be neglected. 

In this subsection we consider two flat models (S model and L model) with 
(i7o,Ao) = (1.0,0) and (0.2,0.8), respectively, and an open model (O model) with 
(0.2,0). The present lengths of the boxes are 



for S, L and O models, respectively. The particle number is 32^ in all models, and 
thus 



respectively. 

The distributions of particles in these models were derived using numerical A^- 
body simulations. The particle size (giving the lens strength) is represented by the 
softening radius (= a(t)xs), which is constant. For the particle size we consider 
two lens models for comparison: 

Lens model 1. All particles in the low-density models {Qq = 0.2) have = 
20/i~^kpc, 20% of the particles in the flat model (1.0, 0) have = 20/i~^kpc, and the 
remaining particles have = 500/i~^kpc. Thus practically jparticles with fi^ = 0.2 
(which we call compact lens objects in previous papersO'EHP) play the role of lens 
objects. Their number density is much larger than the galactic density fig ~ 0.02. 
Lens model 2. 10% of the particles in the low-density model {Qq = 0.2) and 2% 
of the particles in the flat model (1.0,0) have = 20/i~^kpc, while the remaining 
particles have = 500/i~^kpc. Thus only galaxies corresponding to Qg = 0.02 play 
significant roles as lens objects, and the remaining particles are regarded as diffuse 
clouds. 

The lens strength of realistic non-galactic clouds is intermediate between these 
two lens models. The question of which model is better may be answered by obser- 
vational studies involving lens phenomena. 

The time evolution of the distribution of particles was derived by performing the 
A^-body simulation in the tree-code provided by Suto.cHP TheJnitial particle distri- 
butions were derived using Bertschinger's software COSMICS& under the condition 
that their perturbations are given as random fields with the spectrum of cold dark 
matter, their power n is 1, and their normalization is specified as the dispersion 
(78 = 0.94 with the Hubble constant h = 0.5 for (1.0,0) and h = 0.7 for other 
models with Qq = 0.2. 



Lo = a(to)/ = 32.5/i"\ 50/i"\ 50/i"^Mpc 



(3-5) 



m(= PbqLq^ 



/N) = 2.90, 2.11, 2.11 X 10 




(3-6) 
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Now let us consider light propagation described by the null geodesic equation 
with the null condition. Here we use T = | ln[a(t)/a(ti)] as a time variable and 
To = {T)t=to- Then we have dS = 2exp[2(T - To)]dT, so that 



where 
and 



cdt = RcR [f2o + {l-f2o- Xo)S + XoS^]-^/'^e^'^dT, (3-7) 

R = Lo/[{l + zi)N^/'^] (3-8) 
CR = 2ic/Ho)/[Ril + zif% (3-9) 



The line-element is 

ds'^/R^ = -CR'^[f2o + {1-^0- Xo)S + XoS^]-^e^'^{l + a(j))dT^ 

+ (1 - a^)e^^dy^/F{yf, (3-10) 

where y° = T, y' = a{ti)x'/R, (p = (Gm/R)^, Rq = Ra^/ai = (1 + zi)R, 

2Gm 3 Qq 



a = 



and 



F=l- -iRoHo/c)\l -f2o- Ao)(y)'. 



The equations to be solved for light rays are 

dy' 



dT 



= CRe^K\ 



(3-11) 
(3-12) 

(3-13) 



^ = - [3Aoe<^-^o) + (1 - I2o - Ao)]e^(^-^o)^YG(r) + a^K^ 



-^CR-'e^ d<j>/dyyG{T)-2-^K^K^ 

+ {RoHo/cf{l -no- Xo)F-^CRe^_ 
■y^k^K' + ^{l + 2a(f))/G{T) 



where = cr ^dy^dT, 7 = a(cR)^ and 

G(T) = f2o + il-no- Ao)e2(^-^°) + Aoe^C^-^o). 
The null condition is 



(3-14) 



(3-15) 



(3-16) 



While in flat models the solution is straightforward, some care must be taken 
in curved models. The potential is given as a solution of the Poisson equation. 
Because the ratio of the second term to the first term on the right-hand side of Eq. 



102 



K. Tomita, P. Premadi and T. T. Nakamura 



(3-3) is {RqHq/ c)'^{y)'^[Ay /\y\] <C 1 for z < zi, the Poisson equation in a box can be 
approximately expressed as 

P^Vcf = ^^GpB\p{x)/pB - 1], (3-17) 

where we have used F{y) ~ 1 for y ~ [the box size]. F{y) can be approximately 
replaced by the central value F{y^) for y » [the box size], where the index c denotes 
the central value in the box. For point sources with p = mJ2n ^{(^^{v ~ Vn)) ('^ 
the particle number) , we have (j) = (f)i + (1)2, where 

0i = -F(yJe-2^5:^^, (3-18) 

„ \y-yn\ 

and 02 represents the contribution from the homogeneous background density. Here 
let us use for y another vectors y expressing the space in a locally flat manner, where 
the two coordinates are related as 

y= f^dy/Fiy), (3-19) 
Jo 

and the lengths between two points in boxes in two coordinates are approximately 
related as 

Ay = Ay/F{y,). (3-20) 
Then 0i can be expressed in terms of y in the usual manner as 

n\y-y 

Corresponding forces are expressed as fi = d(j)i/dy^ = [d(pi/dy^]/F{y^). It should 
be noted that the contribution of d<p/dT is negligibly small, compared with that of 

In flat models the universe is everywhere covered with periodic boxes continu- 
ously connected as in Fig. 6. In open models it cannot be covered in such a way, 
but we can consider only a set of local periodic boxes connected along each light 
ray, as in Fig. 7. In these boxes we can describe the evolution in the distribution 
of particles in terms of local flat coordinates y, because the size of boxes is much 
smaller than the curvature radius. 

For the integration of the above null-geodesic equations, we calculate the po- 
tential at a flnite number of points on the ray which are given at each time step 
(Z\T). Then particles near one of the points on the ray have a stronger influence 
upon the potential than particles far from any points. To avoid this unbalance in 
the calculation of the potential, we take an average of the potential 0i by integrating 
it analytically over the interval between one of the points and the next point. The 
expression for an averaged potential (p was given in a previous paper. Moreover, to 
take into account the finite particle size as galaxies or non-galactic clouds, we modify 
the above potential for point sources using the softening radii yg = a{ti)xs/R- The 



(3-21) 
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modified potential is produced by replacing (y — by {y — ys)^ + (ys)^ in the 
potential for point sources. 



The initial values of are given so as to satisfy Eq. ( 3-16 ). The integration of 
Eqs. ( pT3D and (3-14) with the modified potential was performed using the Adams 
method. As the time step we used AT = [ln6/2]/A'^5 with Ng = 3000 (in most cases) 
-10000. 

Here we consider the comoving volume of the region through which a ray bundle 
with solid angle 9'^ pass at the interval < z < Zi in the flat model (1.0,0). This 
volume is equal to the volume of the cone Vr = le"^ [2cHo{l-l/^/TTzl)f. Its ratio 
to the volume of the periodic box Vb [= (32/i^^Mpc)'^] is 1.2 x 10^^ for Zi = 5 and 
9 = lOarcsec. This small ratio implies that the influence of the periodicity of the 
box on the the statistics of ray bundles is negligibly small, as long as the ray bundles 
are not directed in some special directions with respect to the box. 




Fig. 6. Light rays and periodic boxes in a flat Fig- 7. Light rays and periodic boxes in an 
space. open space. 



3.2. Statistical behavior of optical scalars 

We treat the deformation of ray bundles over the interval from z = to z = 5 
measured by an observer in a periodic box. Here we consider the ray bundles reaching 
the observer (or emitted backwards) in a regular form such that the rays are put 
in the same separation angle 9, and calculate the relative change in the angular 
positions of the rays which increases with redshift in the past direction. From this 
change we find the behavior of optical scalars. 

Basic ray bundles consist of 5 x 5 rays that are put in a square form with the 
same separation angle 9 = 2 — 360 arcsec. Many bundles coming from all directions 
in the sky are considered. Here we take 200 bundles coming from randomly chosen 
directions for each separation angle. In order to express relative angular positions 
of rays, we use two orthogonal vectors, e*^j^^ and e*^2); in the plane perpendicular 

to the first (fiducial) background ray vector (K^)b- Then the angular coordinates 
[X{m,n),Y{m,n)] of 25 rays relative to the first (fiducial) ray with {m,n) = (1,1) 
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at any epoch are defined by 

X{m,n) = Y,[y\m,n) - y\l,l)] e\^^/yB{l,l) + X{l,l), 

i 

y(m, n) = ^[y^(m, n) - y\l, 1)] e\^^/yB{l, 1) + ^(1, 1), (3-22) 

i 

where y^ = y^ + and ys = [J^iivh)'^]- Since all angular intervals of rays at the 
observer's point are the same (= 9), the differentiation of angular coordinates of 
the rays at any epoch with respect to those at the observer's points is given by the 
differences 

Aii{m, n) = [X{m + 1, n) - X{m, n)]/e, 
Ai2{m, n) = [X{m, n + 1) - X{m, n)]/e, 
A2i{m^ n) = \Y[m + 1, n) — Y{m, n)]/9, 
A22{m, n) = [Y{m, n + 1) - Y{m, n)]/9, (3-23) 

where m and n run from 1 to 5. From the matrix Aij{m,n) we derive the op- 
tical scalars in the standard manner, 0) as the convergence {K{m,n)), the shear 
(7j(m,n), i = 1,2), and the amplification {^{m,n)) defined by 

K(m, n) = 1 — tr(m, ra)/2, 7i(?tt., n) = [^22(^71, n) — ^ii(m, n)]/2, 
72(m,n) = -[Ai2{m,n) + ^2i(m., n)]/2, 

7^ = (71)^ + (72)^ = [tr(m,n)]2 - det(Aij(m, n)), 
/x(m,n) = l/det(^ij(m,n)), (3-24) 

where the trace is tr(m,n) = Aii{m^n) +^22("i)'^)- The average optical quantities 
in each bundle are defined as the averages of optical quantities for all rays in the 
bundle as 



^^K(m,n) /4^ K? = ^Y^{n{m,n)) 



742, (3-25) 



and so on. In this averaging process the contributions from smaller scales can be 
cancelled and smoothed-out. The above optical scalars at the separation angle 6 are 
accordingly derived in the coarse-graining on this smoothing scale. 

For the present statistical analysis we excluded caustic cases, considering only 
the case of weak lensing in the sense of no caustics. The averaging for all non- 
caustic ray bundles is denoted using () as (k^), (7^) and ((^u — 1)^). Because 
K,{m,n), ^i{m,n) and fi{m,n) — 1 take positive and negative values with almost 
equal frequency, (k), (7^) and (fi) — 1 are small. 

In Figs. 8 and 9, we show the behavior of (7^) for various separation angles 
9 = 2 arcsec — 360 arcsec (= 6 arcmin) at epochs z = 1.0 and 2.0, respectively. 
Similarly in Figs. 10 and 11, we show the behavior of ((1 — /i)^) at epochs z = 1.0 
and 2.0, respectively. The calculations were performed at ^ = 2, 5, 10, 30, 60, and 
360 arcsec. The behavior of (k^) is similar to that of (7^), as was shown in previous 
papers. In all figures, the values in the two lens models are shown using solid and 
dotted lines. It is found in the low-density models that the ratios of (k^)^/^, (7^)^^^ 
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and ((1 — /i)^)^^^ in the lens model 1 to those in the lens model 2 are about 4 at 
the separation angle = 2 arcsec. The ratios decrease with the increase of 9, and 
are < 2 at ^ = 360 arcsec. That is, the remarkable difference between the two lens 
models appears at the small angles ^ ~ 2 arcsec. The values of optical scalars for 
9 ~ 360 arcsec^e roughly consistent with the results of Bernardeau et al. (cf. their 
Figs. 3 and 4)eIP and NakamuraE^-* in their treatments. 

At all separation angles and for both lens models, {k?), (7^) and ((1 — /i)^) in the 
open model (Qq = 0.2) are larger than those in the flat model (L). In model S they 
are largest at all angles for the lens model 2 and at large angles {9 ~ 360 arcsec) in 
the lens model 1, but they are smallest at small angles {0 < 30 arcsec) because the 
role of non-galactic clouds on small-scale is small. 

In Tables I, II and III, we show for models S, L and O the numerical values of 
(k), (^2)1/2^ (^)^ ((^_ 1)2)1/2^ ^^^^ j^j^d ^^2^1/2 at epochs z = 1,2,.., 5 for = 2 

arcsec. 

The influence of lensing on source magnitudes is 

Am = ^ log[l + ((1 - ^)2)i/2] ^ i.09((l - ^)2)i/2. (3.26) 

If 6* = 2 arcsec, the value of Am is (0.058, 0.034), (0.077, 0.019) and (0.133,0.022) 
for the lens model (1,2), respectively, at z = 1 in models S, L and O. If the lens 
model 1 is more realistic, therefore, this lens correction may introduce some sensitive 
modification to the standard selection of cosmological models in the [m, z] relation. 
Here it should be noted that the lensing effect at small angles, which is caused by 
small-scale nonlinear inhomogeneities, is important for the magnitude correctiDn of 
small high-redshift objects, such as SNe la, in cosmological observations. ElP "E2P 







Table I. 


Optical quantities in 


model S at S 


= 2 arcsec. 




lens 


z 


{^} 








(7i> 






1 


0.0005 


0.0245 


1.0036 


0.0529 


-0.0032 


0.0251 




2 


0.0004 


0.0398 


1.0075 


0.0858 


-0.0061 


0.0409 


1 


3 


0.0021 


0.0495 


1.0149 


0.1084 


-0.0080 


0.0514 




4 


0.0034 


0.0563 


1.0208 


0.1251 


-0.0096 


0.0593 




5 


0.0042 


0.0617 


1.0253 


0.1390 


-0.0109 


0.0653 




1 


0.0005 


0.0155 


1.0020 


0.0315 


-0.0026 


0.0152 




2 


0.0020 


0.0251 


1.0065 


0.0516 


-0.0049 


0.0238 


2 


3 


0.0030 


0.0305 


1.0096 


0.0636 


-0.0053 


0.0283 




4 


0.0037 


0.0343 


1.0120 


0.0721 


-0.0055 


0.0314 




5 


0.0041 


0.0370 


1.0136 


0.0784 


-0.0055 


0.0339 



For the statistical analysis we used 200 ray bundles reaching an observer in a 
single inhomogeneous model universe. This number of ray bundles may be too small 
to cover the influences from complicated inhomogeneities in all directions. To obtain 
more robust statistical results it may be necessary to use more ray bundles and more 
model universes produced with random numbers. 

Finally, we touch upon recent observations of^osmological shear due to weak 
lensing and their relation to our results. Fort et al.H attempted measurements of a 
coherent shear from foreground mass condensations in the fields of several luminous 
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100 <7S'' 




Fig. 8. The angular dependence of (7^)^''^. Solid and dotted lines represent behavior for lens models 
1 and 2, respectively, at 2: = 1. The separation angle 9 is in units of arcsec. 



Table II. Optical quantities in model L at = 2 arcsec. 



lens 


z 




{^)^/-^ 




((l-M^V^ 


(71 > 






1 


-0.0011 


0.0324 


1.0025 


0.0705 


0.0008 


0.0347 




2 


-0.0036 


0.0675 


1.0134 


0.1572 


0.0018 


0.0714 


1 


3 


-0.0045 


0.0974 


1.0345 


0.2407 


0.0020 


0.1002 




4 


-0.0065 


0.1216 


1.0560 


0.3209 


0.0023 


0.1214 




5 


-0.0070 


0.1405 


1.0805 


0.4000 


0.0030 


0.1372 




1 


0.0001 


0.0084 


1.0004 


0.0174 


-0.0004 


0.0088 




2 


0.0004 


0.0175 


1.0021 


0.0367 


-0.0018 


0.0180 


2 


3 


0.0001 


0.0239 


1.0025 


0.0505 


-0.0020 


0.0245 




4 


0.0001 


0.0291 


1.0038 


0.0618 


-0.0019 


0.0296 




5 


0.0005 


0.0333 


1.0055 


0.0710 


-0.0017 


0.0334 



radio sources. Schneider et al.E^P determined the shear in the field (2 min x2 min) 
containing a radio source PKS1508-05 with z = 1.2, and their result is that the shear 
is about 0.03 for an angular scale of 1 min. This value may be consistent with the 
low-density models in the lens model 1, but more observational data are necessary 
to deduce any conclusions. In order to clarify the lensing strength of non-galactic 
clouds it is important to have observational values of the shear at small angles ^ ~ 2 
arcsec. 
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100 <y^>'' 




Fig. 9. The angular dependence of (7^)^''^. Solid and dotted lines denote behavior for lens models 
1 and 2, respectively, at 2: = 2. The separation angle 9 is in the unit of arcsec. 



Table III. Optical quantities in model O at S = 2 arcsec. 



lens 


z 


{^) 






{(i-^nv- 


(71 > 






1 


-0.0002 


0.0425 


1.0221 


0.1219 


-0.0028 


0.0426 




2 


-0.0054 


0.0796 


0.9961 


0.2210 


-0.0041 


0.0832 


1 


3 


-0.0069 


0.1157 


1.0368 


0.6234 


0.0057 


0.1225 




4 


-0.0081 


0.1519 


1.0790 


0.5638 


0.0078 


0.1602 




5 


-0.0064 


0.1872 


1.1627 


1.1880 


0.0120 


0.1980 




1 


-0.0003 


0.0099 


0.9997 


0.0205 


-0.0005 


0.0120 




2 


-0.0012 


0.0209 


0.9993 


0.0441 


-0.0015 


0.0209 


2 


3 


0.0000 


0.0314 


1.0041 


0.0682 


-0.0035 


0.0284 




4 


0.0010 


0.0410 


1.0089 


0.0915 


-0.0047 


0.0370 




5 


0.0021 


0.0496 


1.0144 


0.1140 


-0.0059 


0.0451 



§4. Analytic evaluation of cosmic shear through the power spectrum 

approach cJ 

In this section we analytically estimate the shear and the Jmageamplification 
caused by large scale density inhomogeneity in the universe. "E^P emanating 
from us to the past, and let x be the affine parameter along the central ray in the 



*' In this section we use units in which c — Hq 
sign convention of curvature tensors. 



— 1 and follow Misner et al.tSi' as regards the 
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Fig. 10. The angular dependence of ((1 — /i)^)^'^^. Solid and dotted lines represent behavior for 
lens models 1 and 2, respectively, at z = 1. The separation angle 9 is in units of arcsec. 



bundle and k'^ = (d/dx)"" be its tangent vector. Deformation in the cross section 
of the bundle is described by optical scalars, the expansion = |Vi/c" and shear 
a = ^e^e^Vi/cb, where = ef + is the complex dyad basis on the cross section. 
The rotation lv = |lm(e*"e^)ViA;ft identically vanishes for a bundle which emanates 
from a single point, so we set uj = in the following. These evolve along the ray 
according to the equations 

^ + 02 + |(t|2 = -7^, & + 2ea = -J^, (4-1) 

where the overdot denotes d/dx, and the sources of deformation are 

and T = ^Cacbd^°'^^k'^k'^- The observable effect of this deformation is the evolution 

of a deviation vector connecting two nearby rays on the cross section: 

i = 9C + aC. (4-2) 

where ^ = CaC"- This deviation vector is written as a linear mapping from the initial 
value of ^ at X = 0: 

e = yl*^o + rto- (4-3) 

Here .^o represents the angular separation between the nearby rays in the observer's 
sky, and Reyl, Imyl and P represent respectively the degrees of focusing, rotation 
and shear in the shape of the cross section as seen from the observer. From equation 
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100 (1 - /if>'' 




Fig. 11. The angular dependence of ((1 — fi)^)^^^. Solid and dotted lines represent behavior for 
lens models 1 and 2, respectively, at z = 2. The separation angle 6 is in the unit of arcsec. 

( pip it follows that 

A = 0A + a*r, A = -{nA + T*r), (4-4) 
r = aA + er, r = -{j^A + nr). (4-5) 

The angular diameter distance D is defined to be (cross sectional area/solid angle 
in the sky)^^ so = (|ylp - |r|2)V2 and satisfies 

D = eD, D = -{n + \a\'^)D. (4-6) 

Initial conditions for these differential equations at x = are = A = r = r = D = 
a and A = D = 1. 

The following form of the metric is used to describe the globally homogeneous 
but locally inhomogenous universe 

ds^ = a^ds^ = a^{rj) [-(1 + 2<P)drj^ + (1 - 2^)-/ijdx'dx^]. (4-7) 

Here a is the scale factor and ^ij is the 3-metric on the constant time hypersurface 
of curvature K = + — 1. The gravitational potential ^ (<C 1) incorporates the 
local inhomogeneity generated by the density contrast 5 = (p — p)/p and satisfies the 
Poisson equation 

V^cp = ^f2,i, (4-8) 
2 a 
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where is the Laplacian associated with (we only consider a sub-horizon scale 
inhomogeneity whose wavelength is much smaller than Since the null 

geodesies are unaffected by the conformal transformation, we work in the scaled 
metric ds^ in what follows. Accordingly the "distances" A and P (Eq.( [4-3D ) in the 
metric ds^ are related to those in ds^ as ^1 = aA and P = aA. The affine parameter 
X is related to the redshift 2: = l/a — las 

X{z)= dx{Qox- Kx^ + \qx'^)-^I'^. (4-9) 



a 



From Eq. (|4J) we obtain TZ = K + V^<P and = e^e-'X^X^-^. We regard 
5JZ = V^'P and &F = eV-'X^Vjf? as perturbations to 1Z = K and ^ = 0, and solve Eqs. 



(4-4) and (|4-5|) peturbatively. The zeroth order solution in an exactly homogenous 



universe is 

K-V2sin(i^i/2^), {K>0) 

A{x) = D{x) = { X, {K = 0) (4-10) 

(-K)-i/2 sinh[(-K)V2;^], < 0) 

and P{x) = 0) ^iid hence the image remains undistorted with only an isotropic 
focusing. The differential equation for the first order quantities 5A = A — A and 
6r = r — P are 

6A = -{K5A + D3JZ), SP = -{K dP + D SF), (4-11) 
which are integrated with the initial conditions Q = 6A = dP = 5A = 6P as, 

5A{x) = - r dx' SJZ[X) D{x') P){X - X'), (4-12) 

JO 

5P{x) = - r dx' &F{x') D{x') D{x - x')- (4-13) 
Jo 

Note that 6A is real, so that there is no rotation in the image to this order. Let us 
define 

K = -5A/D , 7 = 5P/D. (4-14) 

The amplification of the brightness of images relative to that in a homogenous uni- 
verse is given by {D /D)'^ = 1 + 2k to first order, so n (convergence) measures the 
fluctuation in the image brightness due to the lensing by the density inhomogeneity. 
Similiarly 7 (shear) is a dimensionless measure of the image distortion due to the 
lensing. 

The sources 8JZ and 6J- of the image distortion are stochastic variables which 
average to zero: {&Z) = {SF) = 0. To evaluate the statistical properties of this 
distortion, it is convenient to Fourier-transform the density field at time rj as 5{k,r]) = 
J d^x 6{x, rf) exp(— ifc-a;) (valid for k"^ ^ |i^|) and define the power spectrum P{k, rj) 
as {6{k,r])6{k',r])) = {27rfP{k,r])6u{k - k'). Then from Eq. (|]|) the two point 
correlation of 51Z at two locations x and x' is 

{5JZ{x) SJZ{x')) = \qI J ^ eMik ■ {x - x')]. (4-15) 
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Below we calculate the two point correlations of convergence (ka Kb) and shear 
(7a 7b) between images A and B at the same redshifts separated by a small angle 
vector 6 on the sky. Since the integrals in Eqs. ( [4-12D and ( 4-13D can be performed 
along the unperturbed path to first order, k ■ (x — x') in the exponent is written 
as fc^ • 6D{x) + (x ~ x') where fcx (or fcy) is the wavenumber perpendicular (or 
parallel) to the line of sight. When ^ <C 1 we can approximate A^y ^ k±_ and obtain 



{&Z[x)SJZ{x')) = ^-^Ql5^{x 



X) -r 



dk A'^{k,r]) 



k ka'^{i]) 



MkeDix)], (4-16) 



where Z\^(A;, ry) = k^ P{k,r])/{2'ir'^). A similar calculation shows that {SF{x) SF* {x')) = 
{8R.{x)mx')) and {SJZ{x) SF{x')) = 0. From Eqs. (|442| ), ( |443| ) and dtiel) there 
results I 



C(^) = (acaKb) = (7a7b) 



where a(x') is given as Eq. 



^ , . ^=^(xO^^(x-xO 



k ka^ix') 



(4-17) 




Fig. 12. Rms convergence and shear C^^'^{Q) = {k^Y^'^ = (lll^)^''^ as a function of the source 
redshift. The solid, dashed and dotted curves correspond to (i7o, Ao, erg) = (1, 0, 0.5), (0.3, 0, 1), 
(0.3,0.7,1). 

In Fig. 12 the root-mean-square convergence and shear C^/^(0) = {k?)^/"^ = 
(|7p)^/^ is plotted versus the source redshift, and in Fig. 13 the square root of the 
two point correlation C^^'^{6) = (ka^b)^^^ = (7a 7b)^^^ of convergence and shear 
between images A and B at redshift z = 2 is plotted versus the separation angle 
6 on the sky. We use the non-linear power spectrum with the shape parameter 
Q^h = 0.25,Ej^ which reproduces well the nearby galaxy survey data. The solid. 



*•* In conventional units, Eq. (4-17) should be divided by [c/HqY. 
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^ 0.025 - 



separation angle /arcsec 



Fig. 13. Two point correlation C^^'^iO) — (ka k,-b)^^^ = (Ta 7b)^''^ of convergence and shear between 
images A and B at redsliift z = 2 as a function of the separation angle 8 in the sky. The solid, 
dashed and dotted curves correspond to {SIq, Aq, erg) — (1,0, 0.5), (0.3, 0, 1), (0.3, 0.7, 1). 



dashed and dotted curves in the figures correspond respectively to the cosmological 
models with {f2o, Xq, as) = (1,0,0.5), (0.3,0,1), (0.3,0.7,1) (standard, open and 
cosmological-constant models), where the normalizations erg of the power spectrum 
on the 8Mpc//i scale are chosen such that the observed abundance of low-redshift 
clusters is consistent with the prediction of Press-Schechter theory. ® The figures 
represent the result for an infinitesimal source size because we do not smooth the 
density field; for a finite source the density field should be smoothed on the scale of 
the source size. 

Since we are using a power spectrum well-constrained by the observational data 
at 2: ~ 0, the differences between the curves in the figures purely reflect the different 
behavior of backward light propagation into the past in these cosmological models. 
Intuitively, the rms amplification and shear should be larger in higher density uni- 
verses, as explained by the overall factor fi^ in Eq. (4-17). However Fig. 12 shows 
that this effect is overwhelmed by the evolutionary effect of the density inhomogene- 
ity: the growth of the density contrast in the standard (or open) model is fastest (or 
slowest) among the three models according to the gravitational instability picture. 1^ 
Since a fast growth implies that the inhomogeneity is rapidly erased as one goes to 
the past, the lensing effect in the standard (or open) model is smallest (or largest) 
in Fig. 12. Therefore the detection of cosmic shear should allow us to probe the 
growth of the large-scale density inhomogeneity, when combined with observational 
constraints for z ~ 0. On the other hand, Fig. 13 shows that in the open (or stan- 
dard) model the decrease of correlation with increasing angle is fastest (or slowest). 
This is due to the cosmological geometry (Eqs.( p^ ) and ( 4T0| )): for fixed 9 and z 
the physical separation 6D in the open (or standard) model is largest (or smallest) 
among the three models. Thus the slope of the correlation function of cosmic shear, 
dC/dlnO, should be a measure of the geometry of our universe. 
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§5. Concluding remarks 

In §2 it was shown that the multi-lens-plane method can treat many ray bundles 
and realistic galaxy lens models in very useful and practical ways to derive statistical 
behavior of strong and weak lensing owing to the simplicity of the calculations. 
The discussion in §3 shows that the recent progress in supercomputers has made it 
possible to use the direct integration methods for lensing analyses, and that non- 
galactic lens objects on galactic scales (if any) may have more important influences 
on the cosmological lensing (than galactic lenses) on small-angle scales (< 1 min) 
depending on the lensing models. In §4 it was shown that perturbative methods 
based on the CDM power spectrum are very useful for analyses of weak lensing on 
large-angle scales (^ 1 min). 

For the lensing correction to the luminosity of small objects such as super novae 
we must treat ray bundles with small separation angles (< 2 arcsec) for which the 
deflection due to the small-scale matter distribution on galactic scales is important. 
^From the analysis in §3 it is found that the correction can be larger by a factor of 
~ 5 than that due to the perturbative methods. 
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